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1 Introduction 

This paper investigates the role of topological methods in the analysis of canard- 
type periodic and chaotic trajectories. In Sections [1] - [S] we apply topological 
degree [J [2] to the analysis of multi-dimensional canards. This part of the paper 
was written mainly by the first and the last authors. Sections [5] -[3 are devoted 
to an application of a special corollary of the Poincare-Bendixson theorem to 
the existence of periodic two-dimensional canards. This fragment was written 
mainly by the first and the second authors. 

If 14^: M'' ^ R"* is a continuous mapping, C K'' is a bounded open set, 
and y e R"^ does not belong to the image W{dft) of the boundary dfl of fi, 
then the symbol deg(W^, il, y) denotes the topological degree [I] of at y with 
respect to n. If ^ W{dfl), then the integer number "f(W,n) = deg(W, ri,0), 
called the rotation of the vector field W at dVl, is well defined. A detailed 
description of properties of the number j(W, ft) can be found, for example, in 
|2J. In particular, if I denotes the identity mapping, I{x) = x, then the number 
7(J — W, n) measures the algebraic number of fixed points of the mapping W 
in f2. 

Consider the slow-fast system 

X = X{x, y, e) + X{x, y, z, e), 

ey ^Y{x,y,e) -\-Y{x,y,z,e), (1) 
i = Z{x,y,z,e). 

Here 

xeR^ yeR\ zeR'^, 
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and e > is a small parameter. The terms X{x,y,z,e) and Y(x,y,z,e) are 
supposed to be small with respect to the uniform norm: 



snp\X{x,y,z,e)\,sup\Y{x,y,z,e)\ < 1. 



However, no estimates on derivatives of those functions are assumed. 
The subset 



S= {{x,y,z) G X X R'*: Y{x,y,0) = 0} 



(2) 



of the phase space is called a slow surface of the system ([T]): on this surface 
the derivative y of the fast variable is zero, the small parameter e vanishes, and 
there are no disturbances X and Y. The part of S where 



is called attractive {repulsive, respectively). The subsurface L C S which sepa- 
rates attractive and repulsive parts of S is called a turning subsurface. 

Trajectories which at first pass along, and close to, an attractive part of S 
and then continue for a while along the repulsive part of S are called canards 
or duck-trajectories. 

We apply topological decree to prove existence, and to locate with a given 
accuracy periodic and chaotic canards of system ([T]). The canards which may 
be found in this way are topologically robust: they vary only slightly if the 
right hand side of the system is disturbed. This does not mean that the canard 
trajectories are stable in Lyapunov sense. However, unstable periodic canards 
are useful on their own. Whenever an (unstable) periodic canard describes pro- 
cesses which are interesting, for instance, from the technological point of view, 
this process can stabilized using standard feedback control algorithms (note the 
role of the Pyragas control in this area). Topologically robust chaotic canards 
also have a role: their existence implies the existence of an infinite ensemble 
of (unstable) periodic canards. General features of this ensemble and methods 
of accurate localization of each of its members follow from our constructions 
below. Thus, in this case one has a wide choice of possible periodic modes in 
system ([1]), each of which may be stabilized in the usual way. 

2 Periodic canards 

In this section we formulate the main existence result for topologically stable 
canard-type periodic trajectories. 

Assumption 2.1. We suppose that the function X and Z are continuously dif- 
ferentiable, and the function Y is twice continuously differentiable. The func- 
tions X and Y are continuous. 

Emphasize again, that we do not require any smoothness of the functions X 
and Y . In particular, if there is no variable z, then we investigate existence of 





(3) 
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(periodic) canards of the three dimensional system 



X ^ X{x,y,e) + X{x,y,s), 
ey = Y{x,y,e) +Y{x,y,e) 

Here the disturbances X{x,y,s) and Y{x,y,s) are small in the uniform norm, 
but there are no bounds on their derivatives. Even in this three-dimensional sit- 
uation applicability of the standard toolboxes, which are based on asymptotical 
representations of slow integral manifolds, is questionable. 

Loosely speaking, we prove that a periodic canard in disturbed system ([T]) 
exists, providing existence of a periodic canard in the undisturbed system 

x^X{x,y,e), 

[o) 

ey ^ Y[x,y,s). 

A point {xc, yc) is called a critical point of the system ([5]), if it satisfies the 
equations 

(X(Xe,2/e,0),i;'(Xe,yc,0)) =0, (6) 

y(a;„y„0) = 0, (7) 
Yy{x,,y,,Q)^Q. (8) 

This is a system of three equations with three variables, so in general case it is 
expected to have solutions. The existence of critical points is important for the 
phenomenon of canard type solutions, because every canard, which first goes 
along the stable slow integral manifold for t < Q and then along the unstable 
slow integral manifold for t > 0, must pass through a small vicinity of a critical 
point. Thus, we turn our attention to the critical points, and to the behavior 
of ([5]) in a vicinity of such a point. 

A critical point is called non-degenerate, if the following inequalities hold: 

X(Xe,2;e,0)^0, (9) 

y,'(a^c,2/c,0)^0, (10) 
y;;(a;„ye,0)^0. (11) 

We consider only non-degenerate critical points. Note that non-degeneracy is 
stable with respect to small perturbations of the right-hand side of 

To study periodic and chaotic canards of the full system ([1]), we first con- 
sider canards passing through a small vicinity of a critical point {xc,yc) of the 
truncated system ([5|). Without loss of generality we assume that the critical 
point is situated at the origin: 

Xc = yc = 0. (12) 

Consider the auxiliary system 

x^X{x,y,0), 
{{x,y),Y^{x,y))^0. ^ ' 
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If the initial point {xq, yo) lies on the slow surface 

So^{{x,y)eR'':Yix,y,0)^0} (14) 

of the system ([5]) , then (fT3)) is equivalent to 

x = X{x,y,0), 
{x,y) e So. 

The system dTSl) is important because it describes the singular limits of the so- 
lutions of (O which lie on the slow surface. Equations (fT3|) can also be rewritten 
in the following form: 

i = X{x,y,0), 
yYyix,y,0) = -{Xix,y,0),Y,ix,y,0)). 

Due to ([S|) this system has a singularity at the origin. Therefore, the exis- 
tence and uniqueness of a solution of which starts at the origin requires an 
additional assumption and will be discussed in detail later. 

To describe the dynamics of the system ([S]) near the origin, we introduce a 
special coordinate system (a;^^^ , a;^^^ , y) in the three-dimensional space of pairs 
{x,y). We choose x^^^ to be co-directed with the gradient F' ((0, 0), 0, 0), and 
x^"^^ to be orthogonal to x^^^ and y. In the coordinate system {x^^\x^^\y) the 
gradient of Y{x^^'> , x'^^^ , y, 0) at the origin takes the form 

r'(o,o,o,o) = (e,o,o), e>o, (17) 

and the equation ^ takes the form 

=X«(a;(i),a;(2),j/,0), 

=X(2)(a;(i),a;(2),y,0), (18) 
2) = y(a;«,x(2),y,0). 

Equations (flT)) and ([6]) imply 

X(i)(0,0,0,0) = 0. (19) 

Taking into account the non-degeneracy of the origin, we can guarantee the 
inequalities 

x(2)(o,o,o,o) >o, i;';(o,o,o,o) = c>o, (20) 

by changing, if necessary, the directions of the x^^-* and y axes. 

The existence of canards and uniqueness of solutions of pB]) is guaranteed 
by the following assumption. 

Assumption 2.2. 

2x% (0, 0, 0, o)i;';(o, 0, 0, 0) - (o, o, o, o)y;(„^(o, o, o, o) < o, (21) 

X(i)(0,0,0,0) > 0. (22) 
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Lemma 2.1. There exist Ta < < Tr such that system (|13p with the initial 
conditio 

x(0) = y(0) = 0, 

has the unique solution 

w*{t)^{x*{t),y*{t)), Ta<t<Tr, 

and the inequalities 

Yy{x*it),y*it))>0, 0<t<Tr, (23) 
Yy{x*it),y*it)) <0, Ta<t<0 (24) 

hold. In other words, the half -trajectory {x* (t) , y* (t)) , Ta < t < lives on the 
attractive part of the slow surface (jl4p . and the half-trajectory [x* (t) , y* (t)) , 
< t < Tr lives on the repulsive part. 

Lemma 12.11 implies strict limitation on the possible location ol canards of 
system (O that passing near the origin: such canards should follow closely 
the solution w*{t) for a certain interval ta < < tr. The above argument 
shows that a periodic canard should have a segment of fast motion from a 
small neighborhood of some point of the repulsive part of w*(t) to a small 
neighborhood of the attractive part of w*(t); this fast motion is, consequently, 
almost vertical (i.e., almost parallel to the y axis). More precisely, if there is a 
limit of periodic canards as e — > 0, then the limiting closed curve has necessarily 
a vertical segment connecting the repulsive and attractive parts of w* (t) . The 
next assumption ensures a possibility of such vertical jumps. 

Assumption 2.3. The two-dimensional curves Ta, Tr defined by 
Ta = {x*{t):Ta <t<0}, Tr^ {x* (t) : 0<t<Tr} 
intersect, that is, there exist r and a such that 

x*{t) =x*{a)^x*, 

with 

Ta <T <0 < a <Tr. 

Let, for example. 

Then we also require that 

Y{x\y)<Q, y*{T) <y<y*{a). 

To avoid cumbersome derivations, we additionaly require that the curves Ta and 
Tr do not self-intersect. 

Assumption 2.4. The intersection is transversal, that is the vectors x*{t) and 
x*{a) are linearly independent. 
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Now consider the equation 

z = Z*{t,z) = Zix*{t),y*{t),z,0), (25) 

and denote by St the shift operator along the solutions of (|25p by the time T. 

The assumptions listed above are (probably?) known and they guarantee 
existence of (periodic) canards for the system 

Theorem 2.1. Let D he an open hounded set, and let 

Then there exist Eq > and A > such that for any e < Eq and any X, Y 
satisfying 

sup \X{x, y, z, e)|, sup \Y{x, y,z,e)\ < A (26) 

there exists aperiodic solution of ([1]) that passes through the set (Baix*), Ba{y*), D), 
with a going to zero as eo,d go to zero. The minimal period Tmin of this solution 
approaches a — t as Eq, S ^ 0. 

3 Example 

Consider the system 

= -aa;(2) +y/3, 

i;(2) ^ a;(l) + 1^ 

ey = x^-^^ + + x^^^y. 
The curves and intersect transversally on the plane (a;'^\ a;^^^), see Fig. 

m 

0.6- 



0.4 




Figure 1: Curves Fa (solid) and F^ (dashed) for a = 3. 

This system has a periodic canard. Figure [1] graphs the numerical approx- 
imation of this canard, together with the limiting curve, which consists of F^, 
Fr, and a vertical segment connecting them. 
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Figure 2: Periodic canard for a — 3 with e — 0.1 (thick) and the limiting curve. 
According to Theorem 12. II the perturbed system 

±{2) = .x(i) + 1 + 1(2) (x,2;,e), 
ey — x'^^ + y^ + x^^'>y + F(a;, y, s). 

has for a small e a periodic canard for any X^-^\x, y, e), X^^) ^x, y, e), Y(x, y, e), 
which are sufficiently small in the uniform norm. 

4 Chaotic canards 

In this section we study the chaotic behavior of canard- type trajectories of ([!]). 
The method that is used to prove chaoticity combines the scheme suggested by 
P. Zghczyhski [17] with the method of topological shadowing [18| and uses the 
results obtained in [19j . We specifically note that our results require no computer 
assisted proofs, in contrast to typical application of the aforementioned scheme, 
see [ini m] • For another approach in studying of chaos in singularly perturbed 
systems, see |16) and bibliography therein. 

4.1 Definition of chaos 

Important attributes of chaotic behavior of a mapping f-.'R'^ — > R'^ include 
sensitive dependence on initial conditions, an abundance of periodic trajectories 
and an irregular mixing effect describable informally by the existence of a finite 
number of disjoint sets which can be visited by trajectories of / in any prescribed 
order. 

Let U = {Ui, . . . , Um}, m > 1, be a family of disjoint subsets of R*^ and let 
us denote the set of one-sided sequences cu = ujo,uji,. . . by fi^j. Sequences in 
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fl^ will be used to prescribe the order in which sets Ui are to be visited. For 
X £ [JiLi '^6 define I{x) to be the number i satisfying x G Ui. 

Definition 4.1. A mapping / is called U-chaotic, if there exists a compact 
/-invariant set S c[J^Ui with the following properties: 

(pi) for any uj e il^ there exists x & S such that p{x) G C/^^. for i > 1; 

(p2) for any p-periodic sequence uj G there exists a p-periodic point x G S* 
with P{x)eU^r, 

(p3) for each 77 > there exists an uncountable subset 5(77) of S, such that the 
simultaneous relationships 

limsup \l{r{x)) - I{ny))\ > 1, liminf \r{x) - r{y)\ < r, 

i — ^00 * 

hold for all x,y £ S{r]), x ^ y. 

This definition is a special case of iU, A:)-chaoticity from [12] with fc = 1. 

The above defining properties of chaotic behaviour are similar to those in the 
Smale transverse homoclinic trajectory theorem with an important difference 
being that we do not require the existence of an invariant Cantor set. Instead, 
the definition includes property (p2), which is usually a corollary of uniqueness, 
and (p3), which is a form of sensitivity and irregular mixing as in the Li- 
Yorke definition of chaos, with the subset S{rj) corresponding to the Li-Yorke 
scrambled subset Sq. 

4.2 Chaotic behavior 

Let us now change Assumption 1 2 . 31 to the following stronger assumption, which 
ensures the existence of multiple intersections between the curves Ta and F^. 

Assumption 4.1. Let the trajectories Fq and Fj. intersect K > 2 times, that 
is, there exist Ti and (7^, i = 1, . . . , K , such that 

x*{Ti) = x*{ai) = X*, 

with 

Ta < Ti < < ai < Tr, 

and Ti 7^ Tj for i ^ j- We also require that the curves Ta and F^ do not 
self-intersect, and that 

y(x*,2/)<0, ye[y*{n),y*ia,)], i = l,...,K. 

We also change Assumption 12.41 to the following: 

Assumption 4.2. All the intersections are transversal, i = 1, . . . , K . 
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Theorem 4.1. Let Di C M be open, convex and bounded, and let 

Then there exist disjoint sets Hi 3 x*(Ti), Eq > and A > such that for any 
e < £o o-nd any X, Y satisfying 

s\xp\X{x,y,z,e)\,s\xp\Y{x,y,z,e)\ < A 

the appropriately defined Poincare map V '. |Jj Hi x — > of system ([1]) is 
{Hi X Di}-chaotic. 

Let 5 C IJjlli X Di, i = 1,...,K, be the compact 'P-invariant set from 
Definition 14.11 its existence is guaranteed by Theorem 14.11 Denote by £s the 
topological entropy of the Poincare map V with respect to the compact set S, 
see [H], p. 109. 

Corollary 4.1. Under conditions of Theorem \4-l\ for sufficiently small e the 
following inequality holds: 

£s > log A'. 

This corollary follows from the {11^ x D,i}-chaoticity of P and the definition 
of topological entropy, see Proposition 2.1 from [TO]. 

5 Proofs 

5.1 Proof of Lemma 12.11 

The equation has a singularity at the origin, as demonstrated by 
Thus, we only need to show the existence and uniqueness of the solution at the 
origin. To do this, we eliminate the third equation from (|13p by rewriting it in 
a specially selected curvilinear coordinate system in a vicinity of the origin, and 
by proving the one-sided Lipschitz conditions for the transformed system. 

The curvilinear coordinate system (p, q, h) is introduced in the following way: 
we keep the x*-^-* and y axes from the (x*-^-* , y) coordinate system, so that 
they become the p and r axes correspondingly, and direct the q axis along the 
turning line L on the slow surface 5*0. In the new coordinates system ([S]) takes 
the form 

q = Q{p,q,h,e), (27) 
sh = H{p, q, h, e). 

This representation is valid in a sufficiently small vicinity of the origin, which 
we denote by = {{p, q, h) : \p\, \q\, \h\ < So}. Recall that the turning fine of 
the surface 5*0 is described by equations 

Y{xi,X2,y,0) = 0, 
Y'ixi,X2,y,0) = 0. 
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Due to the choice of g, these equations take the form 

p 0, h^O 

in the (p, q, h) coordinates, implying 

i?(0,g,0,0) = i?,;(0,g,0,0) = (28) 

in a sufficiently small vicinity of the origin. Let us now calculate the tangent to 
the turning line at the origin. Due to (fT7)) and (PH)). Y{xi,X2,y,0) admits the 
following representation at zero: 

yixi,X2,y, 0) = S,xi+Cy'^/2 + (pyx2 + axl + bxl + cxiX2+dxiy + o{\\{xi,X2, y)||^), 

where ^ = (0) > 0, C = VyyiO) > 0, = >;"(2)j,(0). Thus, the equations of 
the turning line can be approximated at zero as 

^xi + + ^yx2 + ax2 + bx\ + cxiX2 + dxiy — 0, 

+ ipx2 + dxi = 0. 

From the first equation we have that xi has an order of square, thus it can 
be eliminated from the second equation. Therefore, the tangent vector at the 
origin is 

(0,C,-^). (29) 

Inequalities ([^0]) imply that this vector, and therefore the q axis, forms a sharp 
angle with the the vector w*(0) = (0, X(2)(0, 0, 0, 0), 0), thus HH) and ^ 
become 

P(0,0,0,0) = 0, Q(0,0,0,0) > 0, i?,"^(0,0,0,0) > 0, (30) 

and the Assumption 12.21 translates into 

P^(0, 0, 0, 0) < 0, P^(0, 0, 0, 0) > 0. (31) 

In the new coordinate system (p, g, h) the The attractive part Sa of Sq satisfies 

the equation 

h^-^-pDa{p,q), p<0, (32) 

and the repulsive part Sr satisfies the equation 

^-pDr{p,q), p<0, (33) 

where Da{p, q), Dr{p, q), are smooth functions satisfying Da{0, 0), -0^(0, 0) > 0. 

Substituting r with ([5^ and into (fT^ . we obtain a system of two 
equations describing the half-trajectory of that lives on the attractive part 
of the slow surface: 

p = Pa{p,q) = P{p,q, -^/-pDa{p,q),0), ^^^^ 

q ^ Qa{p, q) = Q{p, q, ~^/-pDa{p, q), 0), 
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and a set of equations describing the half-trajectory of that hvcs on the 
repulsive part of the slow surface: 



p^Pr{p,q)^P{p,q, ^~pDrip,q),0), 

(35) 

q = Qrip, q) = Qip, q, V-pDr{p, q), 0), 

with the initial condition p{0) — q{0) = being the same for both ([34|) and 
([55]) . The existence of solutions of ([M]) and ([55]) follows from the continuity of 
the right-hand side. To establish the uniqueness of solution of (p4)) in negative 
time, we divide the first equation from (|34p by the second, obtaining 



dp 

= Pa{p,q)/Qa{p,q)- 

Then, taking into account that Qo(0, 0) > 0, we prove the one-sided Lipschitz 
condition for Pa{p, q) in the p variable: 

{Pa{Pl,q) - Pa{P2,q)){Pl -P2) > -La{pi -^2)^, 

where —El < Pi,P2 < 0, |<?| < El, and La > 0. This follows from the elementary 
estimate 

d 

-^Pa{p,q) > L, 

which holds in a small vicinity of zero for an appropriate L < 0. This proves 
the uniqueness of the solution p{q), and uniqueness of p{t) and q{t) follows. 
Uniqueness of the solution of psp in positive time is proved in the same way. 



5.2 Proof of Theorem 

Consider the coordinate system (p, q, h) introduced in a vicinity of zero in 
Lemma 12.11 We extend this system to the whole space 'S? by aligning the q 
axis along its tangent vector at zero (|29p outside of this vicinity, and connecting 
to the curvilinear part in a differentiablc way. Thus, we get a global almost 
linear coordinate change, coinciding with the curvilinear change at the origin. 

From this moment, we will be working with this new coordinate system, so 
([5]) takes the form 

P = P{p, q,h,e), 

q^Q{p,q,h,e), (36) 
eh — H{p, q, h, e), 

and takes the form 

p^ P{p,q,h,0), 

q^Q{p,q,h,0), (37) 
{p,q,h) e 5*0. 
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Because the y axis becomes the h axis in the new coordinates, the relationships 
(I23l) and iMl from Lemma [TT] also hold for iMl): 



Hhiw*{t))>0, 0<t<Tr, 
Hh{w*it))<0, Ta<t<0. 

Hence, the attractive and repulsive slow surfaces allow the following represen- 
tation in a certain vicinities of the curves Fa and Tr correspondingly: 

h = ha{p,q), h = hr{p,q), (38) 

where the functions ha and hr are smooth. Denote this vicinities of Fq and F,. by 
fla and Qr correspondingly. In the vicinity f2o of the origin these parametriza- 
tions take the form ([5^ and ([55)1 . 

Substituting ha and hr into (|37p . we obtain autonomous diffcrcntian equa- 
tions which describe the dynamics of the (p, g)-component of the attractive slow 
solutions of the main system ((36)l for negative t: 



p = Paip,q) ^ Pip,q,haip,q),0), ^^^^ 

q = Qa{p, q) = Q{p, q, ha{p, g), 0), 

and of the repulsive slow solutions for positive t: 

p^ Pr {p, q) ^ P{p,q, hr {p, q),Q), ^^^^ 

q = Qr{p, q) = Q{p,q,K{p, q),0), 

These representations are valid in a vicinity of the curves Fq and F^ correspond- 
ingly. In the vicinity of zero ([39| becomes p4|l and (j40| becomes (|35p . 

Consider now a small vicinity of the intersection point (p*, g,) of the curves 
Fa and F^, which existence is guaranteed by Assumption 12.31 Let {po,qo) be a 
point in this vicinity. Denote by to, Po, 9o) = {Pa{t,to,Po,qo),qa{t,to,Po,qo)) 
the solution of (PSjl with the initial condition p(io) = Po, q{to) = 5o, and by 
Wr{t,to,po,qo) — {Pr{'t),qr{t)) the solution of PO)) with the same initial condi- 
tion p{to) = po, q{to) = go- 
Denote 

A = Pa{p*,q*)Qr{p*,q*) - Pr{p*,q*)Qa{p*,q*)- 

Since the intersection between F^, and F^ is transversal according to Assumption 
12.41 A ^ 0, and the numbers u{p, q) and v{p, q) may be defined in this vicinity 

by 

Wa{u{p,q),Q,p,q) e Tr, Wr{v{p,q),0,p,q) G F^. 

Using the coordinates u{p, q) and v{p, q) we introduce for a sufficiently small 
a > a 'parallelogram' set n(Q;), illustrated on Fig. O 

n(a)-{(p,g): \u{p,q)l\v{p,q)\ < a/2}. (41) 

Also introduce the notation for the two 'sides' of this parallelogram: 

R- = {iP,q)- v{p,q) = +a/2sgnA, \u{p,q)\ < a/2}, 
R+ = {{P,q)- v{p,q) = -a/2sgnA, \u{p,q)\ < a/2}. 
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Figure 3: The set n(a) on the plane (p, q). 

The solutions Wa{t,po, qo) = Wait, T,po, qo) of the system (I39|) with the initial 
condition (p(t), (/(t)) = (po, qo) G n have several important properties described 
below. To formulate them, we will need some auxiliary definitions. 

Define the set 

E^{{p,q): (p = OAg<0)V(p<OAq = 0)}, 

which is the union of the left half of the horizontal p coordinate axis and the 
bottom half of the vertical q coordinate axis. Consider a set T{pq, qq) of all time 
moments when the solution Wa{t,po, qo) intersects the set E: 

T{Po,qo) = {t- Wa{t,po,qo) G E}. 

Definition 5.1. Let i{po,qo) be the moment from the set T{po,qo) closest to 
zero, see Fig. H) 

t = argmin \t\. 

t£T{po,qo) 

Below we sometimes omit the point {po, qo) and write simply t, when the argu- 
ments can be uniquely identified from the context. 

Definition 5.2. If i{po,qo) exists, and pa{i,po,qo) = and qa{i,po,qo) < 0, 
then the solution Wa{t,po, qo) together with the point ipo,qo) are called desta- 
bilizing. If Pa{t,po, qo) < and qa(t,po,qo) — 0, then the solution Wa{t,po,qo) 
and the point (^0790) are called stabilizing. 
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Figure 4: The set T and the moment t. 



This definition emphasizes the fact that the solution Ws{t, xo, yo) of the main 
equation ([1]) with a stabihzing initial condition will stay close to the attractive 
half-plane Pa for some time after t > 0, if £ is sufficiently small; if initial 
condition is destabilizing, then h^{t) will rapidly increase after t > 0. The proof 
of this fact will be the subject of several propositions, all leading to Proposition 

Statement 5.1. If a is sufficiently small, then the moment t is defined for all 
{POtIo) G n and is continuous on H with respect to {po,qo). 

Statement 5.2. If a is sufficiently small, then a point {pQ,qo) G 11 is destabi- 
lizing if A ■ v{po,qo) > 0, and stabilizing if A ■ w(po,'?o) < 0. In particular, _R_ 
is destabilizing and is stabilizing. 

Statement 5.3. Let a be suficiently small. If a point {po,qo) is stabilizing, 
then IIfi{wa{t,po,qo)) < for t <t < 6{a) where 5(a) > maxn{HPo: ^o)} > 
depends only on a. 

Statement 15. II follows from the continuous dependence of Wa on (pq, go) and 
from the fact that w* intersects the line q — transversally. 

To prove Statements 15.21 and 15.31 we can consider the projection d{t) of the 
difference between the trajectories Wa{t,po, qa) and w*(t) onto the normal vector 
for w*{t), t > t: 

d{t) = {pa{t) -p*{t))Qa{p*{t),q*{t)) - {qa{t)-q*{t))Pa{p*{t),q*{t)). 

Then the variation of d{t) satisfies the following initial value problem: 

r{t) = {l^Pa{p*{t),q*{t)) + §-^Qa{p*{t),q*{t)))r{t), r{T) = 1, 

and the following equality holds for small Apo = pa — p^,, Ago = qa ^ q*, 
do = ApoQaip*,q*) - AqaPa{p*,q*): 

d{t) = r{t)do + o(Apo) + o{Aqo). 
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Therefore, if do > 0, and Apo and Ago are small, then d{t) > for t < t. 
Assumption 12 . 41 provides that do > for (po, 90) S and do < for (po, 9o) S 
R+. 

Denote by 

We{t,Po,qo) = {pe{t,PQ,qo),qe{t,po,qo),he{t,po,qo)) (42) 
the solution of the system ([551) with the initial condition 

p{t) ^ Po, q{T) = go, Kt) = ^0 = (^*(t) + h*{<^))/'^, 
with (po,9o) e n(Q;). 

Definition 5.3. We say that a is e-close to 6, if a ^ 6 as e ^ 0. 

In our assumptions, the solution Wi,{t,po, go) first rapidly approaches the at- 
tractive half-plane Sa, and then follows Sa until i w 0. The possible subsequent 
behaviors of the solution are the following: 

1. The solution stays close to the attractive part of the slow surface for t < 6, 
where 6 > 0. 

2. The solution begins to rapidly increase in the h direction around t « 0. 

3. The solution follows the repulsive half-plane Sr until a moment t* . More- 
over, for positive t < i* it follows the curve F^. After the moment t* the 
solution may begin to rapidly increase in the h direction, or it may fall 
down back to the attractive part of the slow surface. 

The propositions below provide a more rigorous description of the qualitative 
description above. 

Proposition 5.1. Let e be sufficiently small, and let {wa{t,po, go), ha(wa{t))) € 
Sa for T < t < T (this means that the corresponding solution of (|37p does not 
cross the turning line L). Then the difference h^{t) — ha{pe{t), qe{t)) becomes 
e-small at the moment ti which is e-close to t, and stays e- small for ti < t < T . 

Proof. Lemma |2. II and the choice of a imply that the inequality 

-C2<Kip,q,h,s)<-Ci<0, Ci,C2>0, 

holds in a vicinity fl of the curve {wa{t,po, go), ha{wa{t))), t <t < T. Suppose 
also without loss of generality that all the functions in the right-hand side of 
are bounded by a constant M > along with all the first derivatives. 
By virtue of Assumption 12.31 H{p, q, h, e) < Hq < in a vicinity of the 
vertical line p = p* ,q = q* , h* {t) + Ah < h < ho, where Ah is sufficiently small. 
Thus, the vertical speed of the solution We{t) goes to infinity as e — > 0, and 
therefore the solution reaches the bottom end of this vicinity in an e-small time 
(linear in e). The changes in p^ and q^ are e-small after this time, thus the 
solution exits through the bottom part of this vicinity, entering fi. 
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Denote ip{t) = h^{t) — ha{pe{t), After the solution entered fi, we have 

for ip>0: 

ifit) = kit) - ha{pe{t),qe{t)) = -^H{p,q,h,e) - ^Pe{t) - ^Ut) 

1 h' 

= -Hip, q, h, e) + -j^Pip, q, h, e) + 777 Qb, q, h, e) 
^ H^ 

1 A'P 1 

< -Hip, q,^ + Kip, q), e) + — = -iK- ifiit) + H'^.£) + — 
e Ci e Ci 

By the theorem on differential inequalities, 

ifit) < m = i'Po - eC)e^ + eC, 

where C depends on Ci and M . Calculate the time moment when ifit) becomes 
equal to 2eC: 

-Cjt 

((^0 - £C)e - = eC, 

e (fQ-eC £ Ah* 
m — < -— In ■ 



Ci eC - Ci eC ' 

where Ah* = /i*(cr) — ^*(t). Thus, (fit) gets e-close to zero from above after an 
£-small time interval [t, ii]. 

We obtain the lower bound for (^(t) in the same way: 

1 Af^ 1 

m > -Hip, q,V + haip, q),e)- — ^-iH'^- fit) + H'^.s)- — 
£ Oi £ Oi 

£ Oi 

where (fit) > 0. Therefore 

'Pit) > iPQ + eC)e^ - £C. 

Similar equations can be written for the case ipit) < with C2 instead of Ci. 
In any case, ipit) remains e-close to zero for ti < t < T . □ 

Proposition 5.2. Consider the set C/i(/3) = {ip,q,h) : h > haip,q) — /?} 
defined in the vicinity D,a of the curve for t < by virtue of psp . and in 
the vicinity Qq of the origin for p < by virtue of (|32p . Also introduce the set 
U2i(3) — Hp, q,h) : h> p — (3} in VIq for p> 0, and let Ua = UiL)U2- Denote 

7(e) = inf{/3: Weit,po,qQ) e UaiP) for any (po,9o) e 11}. 

Then 7(e) is £-small. 
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Proof. Suppose that 7(e) is not e-small. Then there exist 70 > 0, e„ 0, 
{Pmq_n), and such that WeitmPn, In) ^ Uai'^o)- Denote the hmit point of 
the sequence We{tn,Pn, Qn) by wq ^ Ua{jo)- Suppose that wq ^ flo, or wq G flo 
and Pa < 0. Consider the scalar product (w^, gradS'a) at the point wq- 

Ci = {we,gTadSa) = -PH'p - QH'^ - \hhI 

Note that H{wo) > 0, i?^(wo) < 0, and the first and second terms are bounded. 
Therefore we can select a sufficiently small e > such that this scalar prod- 
uct is positive, which means that the velocity vector is directed inside the 
set Uai'jo)- Due to the choice of wq we can select n such that e„ < e, and 
We{tn,Pn,<ln) is Sufficiently close to wq, therefore this scalar product should be 
close to Ci by virtue of continuity of the right hand side of and thus 

positive. However, the velocity vector is not directed inside the set Ua{'~io), 
implying that the scalar product must be non-positive, contradicting positive- 
ness of Ci . 

In the case wo € flo, po > 0, the scaler product (w^, grad9J72(7o)) takes the 
form 

C2 = (u;e,grad9C/2(7o)) = -P+^H. 

In this case H{wq) > 0, and the same reasoning applies. This contradiction 
proves that 7(e) is e-small. □ 

This proposition shows that there exists an invariant set e-close to Sa, that 
the trajectories Wa{t,po, qq) of ([55]) with ipo^qa) G 11 do not intersect. By 
applying the same reasoning to solutions of (|36[) in backward time, we can 
obtain the existence of an invariant set Ur e-close to Sr, that the trajectories 
Wa(t,po,qo) which return to the vicinity of the point {p*,q*,hQ) also do not 
intersect: 

Ur = {{p, q,h):h> ha{p, q) - 7(e)} U {(p, q,h) : p > Q ^ h > p - 7(e)}. 

Proposition 5.3. // (xo,?;o) is stabilizing and e is sufficiently small, then 
Pe{t,po,qo) is e-close to pa{t,pa,qo), and qe{t,po,qo) is e-close to (ja(t,Po,9o) 
for all T < t < S{a). 

Proof. According to Statement [Ol Hii{wa{t,po, qo)) < Hq < for t < t < S{a). 
Proposition 15.11 implies in this case that h^{t) — ha{pe{t), qe{t)) is e-small for 
ti < t < S{a), where ti is e-close to r. Therefore, first we have 

\Pe{tl) - Pa{tl)\ < \Pe{tl) ~ Po\ + \pa{tl) - Po\ < 2M{ti - t), 
\qe{ti) ~ qaiti)\ < 2P(<i -r), 

and the statement of the poposition holds for t < t < ti. 
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Denote = \p^{t) - Pa{t)\ + \qe{t) - qa{t)l t > h. Using §^ and 
we obtain 

DB.\Pe - Pal = |P(pe,ge>a(Pe,g£),e) " ^"(^0, ^a, /ia (Pa, ^a), 0) | +P;^7i(£) 

= (P; - Pl,H'jHl,)\p, -Pa\ + {P',- PLH'^/K)\qe -~qa\+ P^s + P^7i(e) 

< C(|pe - Pal + ke - ^a | +7l(£)), 
-Diike - gal < C{\pe - Pa\ + ke - ^a | + 71(e)), 

where D/j denotes the right derivative, 71(e) is e-small, and C is a sufBciently 
large constant. Therefore 

DRijjit) < Ciip{t) + 72(e), V(ii) < 72(e), 

and by the theorem on differential inequalities, 

V'(0<73(e)e^^(*-*^)- 73(e), 

which is e-sniah on a finite time interval ti <t < 5{a). □ 

Proposition 5.4. Let (xo,?;o) be destabilizing and e be sufficiently small. Con- 
sider a small vicinity ft of the point (0, g, 0) — (pait) , qa{t) , 0) , and let t be the 
time when the solution Wa(t) is inside this vicinity. Then Pe{t,pQ,qo) is e-close 
to Pa{t,po,qo), and qeit,po,qo) is e-close to qa{t,Po,qo) for all t < t < i, and 
We(t) exits the invariant set Ur when it leaves fl. 

Proof. By noting that Hh{wa{t,po, qo)) < Hq < for t < t < t and repeating 
the same steps as in Proposition 15.31 wc obtain that Pe(t,pQ,qo) is e-close to 
Pa(i,Po, go), and qe{t,po,qo) is e-close to qa{t,pa,qo) for aU r < i < i. This 
means that We{t) is inside J7 at the moment t. 

The vicinity SIq should be selected small enough that Q > Qo > in this 
vicinity, and a and should be small enough such that C Oq. I|3ip together 
with q < imply that P > Pq > in fi. Let 

^ = {{P, q, h) : \p\ < 6p, \q-q\ < Sg, \h\ < Sh}, 

with sufficiently small Sp, Sq, 6h, such that Sp < -^Sq, and the bottom side of 
Q with h = —Sh is outside Ua, and the top side h — Sh is outside Ur. 

Consider the time moment t when w^(t) leaves this vicinity, ft can do this 
in three possible ways: 

1. p^(t) = Sp. Then h^{t) > -7(e), thus We{t) ^ Ur. 

2. qe{t) = Sq. Due to the choice of Sp, Pe{t) > Sp, thus this case is not 
possible, and case 1 actually takes place. 

3. h^it) — Sh. In this case We{t) ^ Ur by choice of Sh. 

In all cases, leaves C/,-. □ 
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Propositions 15.31 and 15.41 are both valid for sufficiently small e, that is, for 
£ < eo(POj9o)- Due to continuous dependence of the solutions on initial data, 
eo{po, qo) can be selected in such a way that these propositions will be valid for 
£ < So in some vicinity of the point (poi^o)- The next proposition shows, that 
Case 3 above is only possible when the initial condition {po,qo) £ 11 is e-close 
to the curve Ta- 

Proposition 5.5. For every e define a set of points from H such that Proposi- 
tions [573^5.4\ hold for these points with the selected e: 

^{s) = {ipo,qo) e n: £o(po,go) > s}- 

Denote 

7(e) = sup{|u(po, go)| : (Po, 9o) ^ A(£)}. 

Then 7(e) is e -small. In other words, if a trajectory Wg^it^po, qo) does not fulfill 
Provositions \5.3^5.4\ thenv{po,qo) is e -small. 

Proof. Suppose that 7(e) is not e-small. Then there exist 70, e„ 0, (p„, g„) G 
n, such that {pn,qn) ^ A(en)i and \v(jpmqn)\ > 7o- Consider a limit point 
(p,q) of the sequence [^(P:?)! ^ 7o- The latter inequality means that 

{p, q) is either stabilizing or destabilizing, thus Eq (p, q) is defined and positive. 
Select e„ such that e„ < Sq, and (j>n,qn) is in a vicinity of the point {p,q) 
where Eq is defined. Then, on the one hand, we have £„ < Eq, and on the 
other hand, {pn,qn) ^ A(£„), therefore £„ > eo- This contradiction proves the 
proposition. □ 

To prove the existence of periodic canards, we need to define a mapping 
We', n — !■ R^, and to do that we need an auxiliary time moment Se, which is 
associated with the point when the solution either returns to the vicinity of the 
point (p*, (7*, /iq), or deviates from the repulsive curve before or after reaching 
this set. 

Recall that f2o denotes a vicinity of zero, where the coordinate system is 
curvilinear, and relations (|30p - (|3ip are valid, and and Q.^ are vicinities of 
Va and where the representations [551 are defined. By repeating the proofs of 
Propositions 15 . II and 15 .31 in backward time, we can select a constant C such that 
Pe-Pr, qe - qr, K - K < Cs for any initial condition (pe(o-), ^(cr)) = (po, qo) & 
n(2a), he{a) — ho, while ilo- Denote 

n- = {(p, q,r):h< Kip, q) ~ 2Ce} UflrU (R^ \ Qq). 

The definition of the time moment Se{po,qo) is divided into several possible 
cases: 

1. If he{ti,pQ, qo) ^ Ur{2j{e)) for some ti < a + 3a, then Se{po, qo) = (y + 'i.a. 

2. If We{ti,po,qo) € for some to < ti < a — 3a, where to the moment 
when We exits fio, then Se{po, 9o) = o" — 2a. 



19 



3. If We ^ ior Iq <t < a + 3a, then Se{pQ, qa) = a + 2a. 

4. Otherwise, the solution begins to fall back to the attractive part of Sq at 
the moment ti £ [a — 3a, a + 3a]. In this case let Se(poj'Zo) be the first 
moment t2 after ti such that /ie(i2,Po, 9o) = ^Oi or cr + 2a, whichever the 
smallest, or cr — 2a, whichever the greatest. 

Proposition 5.6. If e is sufficiently small, then the moment of time Se{poi<lo) 
is continuous on 11 with respect to pq and qo. 

Proof. Let {po,qo) G 11 and prove that is continuous at the point {po,qo). 
Consider a point {po,qo) close to {po,qo) in all four cases from the definition of 

Se- 

1. /ie(ti,po, go) ^ Ur{2j{s)) for some ti < a + 3a. The intersection between 

and dUr{2j{e)) is transversal, thus by virtue of continuous dependence 
of the solution on initial values there exists a moment ii close to ti such 
that /ie(ti,po, 9o) ^ Ur{2j{e)). Therefore, Se(po,9o) = cr + 2a. 

2. Wi,{ti,po, qo) G for some to < ti < a — 3a, where to the moment 
when We{t,pQ,qQ) exits flo, then Se{po,qo) — o — 2a. If ti = to, then 
We{to,Po, qo) G ^t the moment to when We{t,po,qo) exits fio due to 
continuous dependence on initial values, and s^ijio, qo) = — 2a. 

If ti > to, then the intersection between Wg and lo^ is trasversal at the 
moment ti, thus there exists ti close to ti such that ^^(fijPoj 9o) G ^7 ^ 
and again Se(po7 (?o) = o" ^ 2q;. 

3. Wg{t,po,qo) ^ il^ for <o < ^ < + 3a. Then We{t,po,qo) ^ for 
to ^ ti < (7 + 3a, and Se(po, go) = o" + 2a. 

4. he(t2,poT qo) = /io- The intersection between and the plane h = ho is 
also tranvsersal, as above, thus there exists a moment t2 close to t2 such 
that he{t2,po,qo) = ho. Therefore, s^ is continuous at (po,9o)- 

In all cases, s^ is continuous. □ 

Proposition 5.7. // (poj^o) «s stabilizing, then Se = cr — 2a /or sufficiently 
small e. If {po, qo) is destabilizing, then s^ = a + 2a. 

Proof. Let {po, qo) be stabilizing. According to Proposition l5.31 We(ti) is e-close 
to Waiti) at the moment when exits f^o- Therefore, Wi,{ti) S rt~ , and Case 
2 from the definition of s^ takes place. 

Let now {po,qo) be destabilizing. According to Proposition l5.4[ there exists 
a moment ti such that ^^(ti) Ur{-y{e)). Thus, cither the solution exits the set 
Ur{2'^{e)) and Case 1 takes place, or it remains in the set Ur{2'y{e)) \ Ur{'l{e)), 
thus never entering , and Case 3 takes place. □ 

Proposition 5.8. Let Case 4 from the definition of hold for {po,qo), that 
is, there exists a time moment t2 close to cr such that /i£(t2,P0i ^o) = ^o- Then 
v{po,qo) and u{pe{t2),qe(t2)) are e-small, and t2 is £-close to a + u{po,qo) — 

V{pe{t2),qe{t2)). 
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Proof. The fact that v{pq, qq) is e-smah is a direct corollary of Proposition 15.51 
and e-smallness of u{pe{t2), 9e(i2)) is an equivalent statement in backward time. 

Consider now the vicinity fig of the origin. The solution enters it e-close 
to Fa, which has 5 < at the point of entry, and exits e-close to r,., which has 
g > 0. In this vicinity Q > 0, thus there exists a unique time moment t when 
Qe — Q,w^e ilo- It is sufficient to show that i is e-close to u{po, qo); e-closeness 
oi t2 — t to a — v{pi;{t2), qe{t2y) can be shown by repeating the same steps for 
the system in backward time. 

Suppose the contrary, that is, there exists a sequence £„ 0, {pn,qn) G 
11(0;), such that |in — u(po, 9o)| > ^3 > 0. Let {p,q,t) be a limit point of 
ipn, qn, tn), and let, for example, t < u{po, qo) — t^. Then, according to Proposi- 
tion l5.3l applied at the moment t, qe{t,p, q) should be e-close to q*{t+u(j)o, qo)) < 
0, however, for sufficiently small Sn and (pn, qn,tn) sufficiently close to {p,q,t), 
qeitmPm qn) — 0, thus we arrive at a contradiction. In the case t > u{po, qo) + t3 
we apply Proposition 15.31 to an appropriate time moment < u{po,qo) where 
q*{t4 + u{po, qo)) > — C, and qe{tA,Pn,qn) < this moment exists because 
< Qi < Q{p,q,h) < Q2 inn. □ 

Using the moment Si.(po,qo), we define a mapping of the set Il{a) into 
the plane {p, q). The definition is divided into the following two cases: 

Case 1. If 

a - < Si,{po,po) < (t + 3a/2, 

then 

W^{PO,Po) = (Pe(Se,?30,9o), PO,9o))- 

If we identify the plane (p, q) with the two-dimensional subspace 

Po = {{p,q, ho):p,qeR} 

of the phase space of system (l36l) . then in Case 1 the value We{po, qo) coincides 
with the intersection of the trajectory (j42|) with Pq, as long as the corresponding 
intersection time is close to a. 
Case 2. If 

3a/2 < \se{po,qo) ~ <j\ < 2a, 

then 

... , , 2q! — Isj — (t| Ise — ctI — 3a/2 

We{po,qo) = 77; [Pe[Se,PO,qo):qe{Se,PO,qo))^ W (Se). 

a/2 a/2 

This means that Wdpo, ?o) is a continuous convex combination of the intersec- 
tion of the trajectory ((42|) with Pq and the point w*(se) as long as the discrep- 
ancy \se — cr\ between the corresponding intersection moment Sg and a in the 
the range from 3a/2 to 2a. Moreover, if the discrepancy equals a, then the 
definition is consistent with Case 1; if the intersection moment equals a ± 2a, 
then We coincides with w*{a ± 2a). 
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Lemma 5.1. If e is sufficiently small, then mapping W^{xo,yo) is continuous 
with respect to {xo,yo). 

Proof. Follows from the definition of and Proposition 15.61 □ 

The next lemma establish correspondence between the fixed points of 
and the periodic solutions of (j36|l . 

Lemma 5.2. Let {p,q) G n(Q;) be a fixed point of the mapping Wg. Then the 
solution We{t,p,q) is periodic. 

Proof. If suffices to show that if (p, q) is a fixed point, then Case 1 from the 
definition of holds, so that the value of corresponds to a point on the 
trajectory We{t,p,q) and is not adjusted, as in Case 2. Note that if Cases 1-3 
from the definition of s^ take place, then We{p,q) ^ H. Thus, there exists a 
moment t2 close to a such that h^{t2,p,q) = ho. The same argument shows 
that Se — t2, otherwise We {p, q) ^ H. 

Suppose that Case 2 holds, and let, for example, Se{p, q) > (T + 3a/2. Denote 
{PjQ) — WeiscP, q). Proposition 1 5 . 81 shows that v(j),q) and u{p,q) are e-small, 
and 

3a/2 < Se(p, q) - a = u{x, y) - v{x, y) + 7(e), 
where 7(e) denotes an £-small value. Recall that \u{x,y)\ < a/2. Thus, 

v{x,y) < —a + 7(2) < —a/2. 

Note that the linear combination in the definition of moves the point (p, q) in 
the direction of the point (0, —2a) in (u, u)-coordinates, thus further decreasing 
v{p,q), therefore we obtain v{p,q) < —a/2, which is impossible because v{p,q) 
must be e-small. This contradiction proves that only Case 1 can hold for {p,q). 

□ 

Finally, we calculate the rotation of the vector field / — We on n(a). 

Lemma 5.3. For sufficiently small e the rotation 7(/ — We, n(a)) of the vector 
field id — We at the boundary of the set 11(a) is defined by 

7(/-W^e,n(a)) = sgn(A). 

Proof. Let for example, 

A>0. 

Consider i?_ and the upper and lower sides of the parallelogram 11(a). 
Propositions 15.41 and 15.31 imply that 

We(p,q) ^w*(<T + 2a), (p,q)£R-, and 
We(p,q)=w*(a~2a), (p,q)eR+. 
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In other words, 

u(W,{p,q)) = 2a, {p,q)eR- (43) 
u{We{p,q)) = -2a, {p,q)eR+. (44) 

Also, Proposition 15.81 implies that 

lim v{We{p, q)) ^0. (45) 

The relationships (|^5)l ~ (P5|) imply that in the coordinates {u, v) the mapping 
We on the boundary of 11 is close to the linear mapping Li{u,v) = (0, — 4u), 
thus the mapping / — is close to (and therefore co-directed with) the linear 
mapping L2{u,v) = (u, 5w), and the result follows from the properties of the 
rotation number. □ 

Lemma 15.31 implies that the mapping has a fixed point [p, q) on the 
set n(a), which defines a periodic solution of (|36p according to Lemma 15.21 
Moreover, according to Proposition 15. 8[ both v{p,q) and u{p,q) are e-small, 
and consequently (p, q) is e-close to {p* ,(?*), and also the minimal period of the 
solution starting from (p, q, ho) is e-close to tr — t. Thus, we have proved the 
following Statement: 

Statement 5.4. For any sufficiently small e > there exists a periodic solution 
of system (|36p . and thus (O, that passes e-close to the point (x*,y*). The 
minimal period T^ihi of this solution approaches a ~ t as e — > 0. 

Consider now the system ([Ij where X, Y satisfy ([^5]) with a sufficiently small 
A > and an initial condition (poi 90: 'ioj ^o) with (po:9o) S n(a), zq G D. By 
repeating the four cases in the definition of the moment Se{pQ, qo), we now de- 
fine the moment Se,\{po,qo, zq) which has the same properties as Se{po,qo)- 
Similarly, we define the mapping We,x{po,qo, zq): H x I) — > R^. For suffi- 
ciently small e, A this mapping is continuous. Introduce also the mapping 
Ve,\{pQ, qo, zq) : UxD R'' which is defined as follows. Let z{t;po, qo, ho, zq) de- 
note the z-componcnt of the solution of the system ([T]) at the corresponding ini- 
tial conditions. We define 14,a(po: Qo, zo) as closest to z{se,\{po, qo, zo);po, qo, ho, zq) 
point which belongs to the ball of the radius Ca centred at Scr-rizo) with an 
appropriate constant C. 

By construction for small a,e,X the vector field 

/-(W^e,A,l^s.A) (46) 

is not anti-directed to the vector field 

at the boundary of the domain 11 x Z?. Moreover for small a,e,X each singular 
point of the field (|46)) generates a required periodic solution. It remains to apply 
the product theorem, by which 

7(7 - (We,A, K,a), n X D) = 7(7 - S^-r, D) sgn(yl) ^ 0. 

Thus, Theorem 12. 1 1 is proved. 
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5.3 Proof of Theorem 

Here we outline the proof of the Theorem 14. II for the case K = 2. It shares many 
parts of the results from the previous subsection, so we will only describe the 
main changes. The proof for arbitrary values of K is obtained by appropriately 
adjusting the definitions of and below. 

First, we introduce the local coordinates Ui, Vi in the vicinities of the points 
x*{Ti), and define the 'rectangular' sets I\-i{a) = {{ui.Vi): \ui\, \vi\ < 012}. The 
coordinates Ui, Vi can be extended onto a vicinity D Ilj of the curve in 
such a way that the transformations hi : (ui,Vi) ^ {p, q) are homeomorphisms. 
In the coordinates Ui,Vi all the sets 11^ have the same representation 11* = 
{{u,v) : \u\, \v\ < a/2}. 

Suppose that cti < a2- Consider a trajectory w^it^pQ, go) of the system (|36p 
with the initial data (po, lo, hg^i) G Ilj, i — 1,2, where ho^i = {h* (Ti) + h* (ai)) /2. 
Define the time moment in the following way: 

1. If h^{ti,pQ,qo) ^ Ur{2-f{e)) for some ti < (T2 + Sa, then sf^{po,qo) = 
(72 + 2q;. 

2. If Welti, po,qo) £ rij for some to < ti < ai — 3a, where to the moment 
when We exits Oq, then {po, go) = ~ 2q:. 

3. If We ^e fo'^ ^0 < ^ < o'2 + 3a, then sf{pQ, go) = <y2 + 2a. 

4. Otherwise, the solution begins to fall back to the attractive part of 5*0 at 
the moment ti £ [ui — 3a, 02 + 3a]. Consider the following two subcases: 

(a) If ti G [uj — 3a, (Tj + 3a], let {po, go) be the first moment ^2 after 
ti such that he{t2,po, qo) — ^o.jj or + 2a, whichever the smallest, 
or Oj — 2a, whichever the greatest. 

(b) Otherwise, let s^{po,qo) = (a2 — g\— 4a) (ii — a\~ 3a)/ ((T2 — <j\ — 
6a) + (7i + 2a. 

This moment in time is continuous with respect to {po,qo)- Let us now 
introduce the Poincare map of the system ([36]) . 

Definition 5.4. If Case 4(a) above holds for sf^ , and he{sf ,pq, go) — hoj, then 
the Poincare map Vpq is defined in the point [pq, go) by 

'Ppq{PQ,qo) = {Pe{Se ,PQ,qo),qe{Se ,PQ,qo))- 

Now define the mapping : IJ^ Ilj ^ M^: 
Case 1. If 

(Ti - 3a/2 < sf (po,Po) < CT2 + 3a/2, 

then 

Wf{PO,Po) = {Pe{Se ,PO,qo),qe{Se ,PO,qo))- 

Case 2. If 

(71 - 2a < sf (po, go) < 0-1 - 3a/2, 
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then 

rTrK / \ 2a — Isf- - (Ti| , , ^ , K XX \sf- ~ (tA — ia 12 

W^[po,qo) = 4- ^(pe(sf ,PO,go),ge(sf ,PO,go)) + ^ JT, ( 



Case 3. If 
then 



Qf/Z a/Z 

(72 + 3a/2 < sf (po, go) < + 2a, 



rr.K/ N 2a-|sf-(T2| ;f ^ , |gf -cr2|-3a/2 

a/2 a/z 

The mapping is continuous with respect to (po 190)1 ^-^d the following 
counterpart of the Lemma 15.21 holds : 



Lemma 5.4. Let (po, qo) £ Ili and We(po, 90) — (Pi <?) £ Hj - T/ien i/ie Poincare 
map Vpq is defined in the point (po, go), and Vpq{po, lo) = ip, <?)• 

To prove chaoticity of Vpq, we use the notion of (F, VF)-hyperbolicity from 

m- 

Fix two positive integers du, dg with d^ + ds = d. Let V and W be bounded, 
open and convex product-sets 

V ^ F^'') X F^^) C M*' X W ^ VK^") X c M'^" X 

satisfying the inclusions £ W, and let g : V M''" x W^' be a continuous 
mapping. It is convenient to treat g as the pair [g^"^^ , g*-"') where g'-"^ : V — > R'^" 
and g^") : K''^ . 

Definition 5.5. The mapping g is {V,W)-hyperbolic, if the equations 



X F^^') fl W^"' = 0, .g(F) fl (W^"' X (M'^= \ M^(^))) - 

(47) 

hold, and 

^(g(«)^y(«))^0. (48) 



Here S denotes the closure of a set S. 



The first relationship (|47p means geometrically that the image of the 'it- 
boundary' dV^"'' X y of docs not intersect the infinite cylinder C = W ^ x 
R''= ; analogously, the second equality means that the image of the whole 

set g(V) can intersect the cylinder C only by its central fragment " x W^^\ 
Thus the first equation (|47p means that the mapping expands in a rather weak 
sense along the first coordinate in the Cartesian product R''" x R'*", whereas 
the second one confers a type of contraction along the second coordinate (the 
indices '(«)' and '(s)' refer to the adjectives 'stable' and 'unstable'). 
The following theorem follows from Corollary 3.1 in [TO] : 
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Theorem 5.1. Let f : — > R'* be a continuous mapping. Let there exist 
homeomorphisms hi and product sets Vi such that hj^fhi is {Vi,Vj) -hyperbolic 
for all i,j, and let the family U of connected components of the union set [J /li(l^) 
have more than one element. Then the mapping f is lA-chaotic. 

The following lemma establishes the (11^, nj)-hyperbolicity of . 

Lemma 5.5. For every sufficiently small e the mappings — hJ^W^hi 

are {li^^IV^) -hyperbolic. 

Lemma 15.51 shows that conditions of the Theorem 15.11 are satisfied, thus the 
mapping is ^Y-chaotic. The set S in our definition of chaos is invariant, thus 
on this set coincides with Vpq by virtue of Lemma [STU This establishes the 
chaotieity of Vpq. 

Finally, we consider the full system ([1]) and define the time moment s^^(po, 90,-^0)) 
the Poincaremap V{po, qo, zq) : and the mapping Wf^{po, qo, zq) ■ IJi xZ)i 
R^. These definitions repeat almost literally those of s^ , Vpq and with ap- 
propriate modifications. 

The conditions of the Theorem require that Sa.-nDi C Dj, which together 
with Lemma 15.51 implv that mapping 

{W^,x,z{se^\{po, qo, zq);pq, qo, zq)), 

which is contracting along the (u,j,z) coordinates and expanding along the Vi 
coordinate, is (11, x I?i,n, x _Dj)-hyperbolic for all i,j. This implies that the 
Poincare map V is chaotic. Thus, Theorem 14. II is proved. 

6 A corollary of the Poincare-Bendixson theo- 
rem and periodic canards 

Consider the differential equation 

X = f (x, a) (49) 

with a smooth f. Here x — {x, y) e R^, and a £ [a_, a+] is a parameter. Let F 
be a Jordanian curve which bounds the open domain D. Suppose that for any 
a £ [a_,a+] there exists a unique equilibrium ea G -DljF, and 

det J(ea) > 0, a- < a < a+, (50) 

where J denotes the Jacobian. We also suppose that 

tr J(e^_) -tr J(e^+) < 0, (51) 

where tr J stands for the trace of the Jacobian. 

Proposition 6.1. Let .system (j49[) has no cycles confined in D[JT for a = 

a-,a^. Then for some a £ (a_, a+) there exists a cycle of system (|49p which is 
confined in D[JT, and which touches F. 
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Of course, the gist of this statement is in the last three words: . . which 
touches r ". This proposition is a corollary of the Poincare-Bendixson theorem, 
see the next section for a proof. 

We present below four simple examples to demonstrate the role of Proposi- 
tion [HIT] in analysis of periodic two-dimensional canards. 



Example 1. Consider the system 

i^y, ey^~x + F{y + a) (52) 

with a small positive e. Suppose that F{0) ~ 0, F'{y) < for y < and 
F'{y) < for y > 0. The curve x — F{y) is a slow curve of system (|52p for 
a = 0. The branch x — F(y), y < 0, is the attractive part of the slow curve, 
and the branch x — F[y), ?/ > 0, is the repulsive part. The origin is the turning 
point. Periodic canards are periodic solutions of system (j52p which follow for a 
substantial distance the repulsive branch, see Figure We say that at a = 
system (|52p has a periodic canard of magnitude a > 0, if to any small e > 
one can correspond and a periodic solution {x^,a^{t),ye.a^{t)) of the system 
X = y, ey = —x + F{y + a^), such that: 

ma.x{xs^aAt) ■ ye,aAi) = 0} = a. (53) 
In our case a periodic solution may visit the upper half-plane y > only 




Attractive branch 



Figure 5: Attractive and repulsive branches of the slow curve, and an example 
of a periodic canard (dotted line). 

traveling along the repulsive branch of the slow curve. Thus, this definition is 
consistent with the informal explanation given above. 
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Proposition 6.2. There exists a periodic canard of system (j52p of any given 
magnitude q > 0. 

Proof. Note that for any value of a the only equilibrium is given by 

e, = (F(a),0). (54) 

Thus 

J(ea) = 
and the inequalities 

det J(ea) 1/e > 0, tr J(ea) = F'(a)/£ < (>)0 for a < (>)0 (55) 



1 

-1/e F'{a)/e 



follow. 

We choose as the boundaries a_ < < any numbers satisfying 
-a + F(a_)<0 and - a + F(a+) < 0. 



(56) 




Figure 6: The domain D is bounded by the curve T. The triangle A(q;,/3) 
belongs to the area where —x + F{y + a) < for all a e [a_, a+]. The domain 
D contains all slow curves x — F{y + a), a; G [0, a G [a_, a+]. 

Let us describe the domain D, see Figure [51 We choose a number (3 > a 
such that the triangle 



A(a, I3) = {{x,y):a<x<l3,\y\<x-a} 
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belongs to the area where + F(y + a) < for a_ < a < a+. In particular, 

-x + F{y + a)<Q for {x,y) e A{a, (3), a e [a_,a+]. (57) 

Existence of such f3 follows from ([55]) . We choose also the numbers C- < < C+ 
satisfying 

-x + F{(±+a)>0 for x < /3, a_ < a < a+. (58) 

Consider the open quadrangle Q which is bounded from south and from north 
by the lines y = C,- and y = C,^, bounded from east by the line x — [3, and 
bounded from south-west by the line x + y = C- • Denote by D the open set 
Q \ A(a, (3), and denote by F the boundary of D. 

Note that the equihbria Gq, a G [a_, a+], belong to D |J F by (|5^ : thus a_ < 
< a+ and (|55|) guarantee that ([50)1 and ([?T|) hold. To apply Proposition lG.ll it 
remains to show that there are no cycles confined in -D IJ F for a € {a-, a+}- 





V 




-0/2 

















Figure 7: The sub-domain 13* C 13 is bounded by the bold dashed line. A 
trajectory which is confined in D cannot leave the domain D^^ once it entered 
D*. Each periodic solution x*(t) = (a;*(i),y*(i)) which is confined in D must 
visit 13*, because it must cross the bold segment of the axis = 0. The areas 
are shrinking within Z3*, and thus there are no cycles there. 

Consider the case a = a_. Introduce the auxiliary sub-domain C D 
which is bounded from north by the line ?/ + a_ =0, from west by the line x = 
F{a-/2) and from south-cast by the graph of the function —x+F{y+a-/2) = 0, 
see Figure [7l For small e a trajectory 'x.^{t) = {x^{t)^y^[t)) which is confined 
in D cannot leave the domain D^, once it entered D*. To prove this claim, we 
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note that for the small e the velocity vectors x point inward at the part of the 
boundary of -D* which belongs to D. Indeed, at the west boundary we have x = 
y > 0; at the north boundary the inequality y — {—x + F{y + a_)/e < holds, 
and at the south-west boundary the velocity vectors point almost vertically up 
for small e. Moreover, each periodic solution x*(t) = ?/*(t)) which is 

confined in D must visit D*. Indeed, because x ^ y, the solution x*(t) must 
visit both half-plane y < and half-plane y > 0. Thus x*(i) must cross some 
times the axis y = from above, i.e., for x > Xa_ ~ F{a^/2)] it remains to 
note that the whole interval 

{(a;,0) : F{a-/2) < x < a} 

belongs to D*. 

By the italicized parts of the previous paragraph, each cycle which is confined 
in D must be confined in D^. However, within _D» the inequality tr J(x,?/) = 
F'{y+aJ)/e < holds, the areas are shrinking, and therefore there are no cycles 
there. The case a = a_ is completed, and the case a = a+ can be considered 
analogously in the backward time. 



t t t 
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F(y+a_)=x 
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1 f(y+aj=x 1 - 





Figure 8: A periodic orbit which is confined in D IJ F may touch F only at the 
point (a,0), because at all other points of F at least one end of the velocity 
vector points strictly outward F. 

Thus, by Proposition l6.1[ for any small e > there exists a periodic solution 
(s^e.oe (Oi J/e.oe (^)) whose trajectory is confined in _D1JF and touches F. On the 
other hand, a periodic orbit which is confined in 13 (J F may touch F only at the 
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point (a,0), see Figure [5] at any other point at least one end of the velocity 
vector points strictly outward F. (At the north, south and south-west parts of 
the boundary this is due to almost upward orientation of x for small e, see ([57)) : 
at the sides of the triangle A(q;, (3) (apart of the point (0, a)) — due to almost 
downward orientation of x, see (|58p ; and at the vertical fragments of the east 
boundary — due to x — y ^ Q.) Thus, the family {Xf^^a^ (t),ye,ae (t)) is a periodic 
canard of the required magnitude a, and the proof is completed. 




Figure 9: Attractive and repulsive branches of the slow curve, and an example 
of a periodic canard (dotted hue). 

Statements similar to Proposition 16. 21 provide no information about asymp- 
totic of Oe, and on stability of canards. Still they could be useful in applications: 
the canards which existence is known can be further located and stabilized via 
a suitable feedback in a usual way. Note also that we do not guarantee that the 
a canard of the magnitude a has only one jump point per period the structure 
of a canard may be trickier, see Figure [HI 

Example 2. Consider system (|5^ with a bimodal function /. Suppose that 
F{0) = 0, F'{y) < for y < and for y > > 0, whereas F'{y) > for 
< y < fi. The curve x = F{y) is a slow curve of system (|52|) for a = 0. 
In particular, the branches x = F{y), y < 0, and x = F{y), y > ^, are the 
attractive parts of the slow curve, and the branch x — F{y), < y < fi, is the 
repulsive part. The origin and the point (F{^),fi) are the turning points. This 
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modification of the first example is similar to the classical Lienard equation. 



y 




F(y)=x 



Figure 10: Early canard (dotted) and late canard (dashed for a bi-modal func- 
tion F{x). 

Periodic canards of a magnitude a may exist only for < a < F{fi). More- 
over, there are two possible structures of a canard: solution may jump down, or 
jump up from the unstable part of the slow curve. We will use wordings early 
and late canards correspondingly. Let us give formal definitions. 

We say that at a = system ([52l) has an early periodic canard of magnitude 
a > 0, if the relationship ([53]) holds, and we say that the system has a late 
periodic canard of magnitude a > if instead 



Proposition 6.3. There exists an early and a late periodic canards of any given 
magnitude a £ (0, -F(/i)). 

Proof is similar to that of Proposition 16.11 As a± one can chose any small 
numbers satisfying a_ < < a+; the inequalities ([50|) . ([51]) are evident. A 
possible construction of the region D is given in Figure lllb . Non-existence of 
cycles at a = a_ may be proven as before. For non-existence of cycles at a = 
see Figure [TTb. 

Example 3. As the next example we consider system (l52t where F{y) is a 
continuous piece- wise monotone function which satisfies the following conditions 



= a. 



(59) 



F(0) = 0, F{y) > 0, y ^ 0, lim F{y) = ±oo. 



(60) 
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F(y)=x 

a) 



F(y+a j=x 

b) 



Figure 11: a) Schematic image of a suitable region D bounded by the curve T. 
b) To prove non-existence of cycles at a = we consider region D^,. This 
region is a repeller, and it doesn't embrace any cycle, since the areas are growing 
within. Thus any cycle should cross the horizontal axis to the right of D^, and 
for small e must be close to the relaxation cycle (dotted) . However for small a+ 
this relaxation cycle does not belong to the region D. 

We also suppose that al local extrema of this function are pairwise different. 

For a given Xq > F{yo) we say that that at a = system (152]) has an early 
(xq, yo)-periodic canard, if to any small e > one can correspond and a 
periodic solution {xg^a^ (t), y^a^ {t)) of the system x = y, ey — —x + F{y + a^), 
such that 

max{a;e,a, (0 : Ve^aA^) = Vo} = xo- (61) 

The late canards for the case xq < F{yo) are defined analogously, with the 
difference that (pT|) is swapped by 

min{a;e,a, (t) : J/e.a. (t) = yo} ^ xq. (62) 
Introduce the auxiliary function 




min^>j,F(z) y > 0, 
max^<j,F(z) y < 0. 



Proposition 6.4. There exists an early {xo^yo) -periodic canard for any xq > 
F(yo, and there exists a late {xQ^yo) -periodic canard for any F{yo) < Xq < 
F*{yo)- 

The proof combines the proofs of two previous propositions. 
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a) 



b) 



Figure 12: a) An example of the "super-late" canard for a multi-mode function 



b) The line x — F*{y) (bold) versus the line x = F{y). Early (xq, ?;o)-canards 
exist for any point (xq, yo) located to the right of the line x = F(i/); late {xq, yo)- 
canards exist for any point (xo^yo) located strictly between the lines x = F{y) 
and X = F*{y). 

Example 4. As the last example we consider the system 



X = F{x, y) = x{p - f{y)), ey = G{x, y, a) ^ y{~q + x{r + g{y) - ah{y))). 

(63) 



Here p,q,r > are given numbers, /(O) = .g(0) = h{0) = 0, f'iy),g'iy), h'{y) > 
for y > 0, e is small, and a is a parameter. Suppose also that 



System (|63p has been recently used in population dynamics. Loosely speaking, 
the functions r -I- g{y) and h{y) describe facilitation and competition between 
predators respectively. The equation (|64p means that the competition prevails 
for denser populations of predators. An instructive example of the functions 
9{y), Kv) is given by 



g{y) = aiy + ^22/' + • • • + amy™, h{y) = (3nr+^ + /32y"+' + • • ■ + /?ny'"+", 

(65) 



where all coefficients are non-negative, and at least on ai and at least one Pj is 
strictly positive. Loosely speaking, Ui measure intensity of mutual facilitation 



ny)- 



lim 9{y)/h{y) = 0. 



(64) 
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between i + 1 predators, whereas Pj measure intensity of competition between 
m + i + 1 predators. Another similar example is given by 

I'M i-N 

aiy) = / v{a) da, h{y) = / w{a) da, (66) 
Jo Jm 

where the weight functions v{a), w{a) are positive and bounded, and < M < 
N. 

The system of equation to find "canard-susceptible" triplets (a;*,?;*, a,) is 

F{x,y,a)=0, G{x,y,a)^Q, G'y{x,y,a) ^ 0. 

In the positive quadrant x,y > this can be rewritten as 

f{y)^a, x{r + g{y) - ah{y)) = q, g'{y)=ah\y). 

Since /(O) — and f'{y) > for y > 0, there exists a unique y* > which 
satisfies /(y) = a; thus a* = g' (y*)/ h' (y^,), and a;* = q/{r + g{y*) — a*h{y^,)). 
We suppose that a;* is positive, that is that the inequality 

r + g{y*) > a^h{y^) 

holds. 

In the positive quadrant the slow curve is given by 

x = X{y)^q/{r + g{y)~a,h{y)), < y < 7], (67) 

To avoid non-principal complications we suppose that the function g'{y)/h'{y), 
strictly decreases for y > 0; this is always true in the case ((65)) or ([66)1 . (For 
instance, in the case (I55t we rewrite g'{y)/h'{y) as gi{y)/hi{y) with gi{y) = 
9'i.y)/y"^,hi{y) = h'{y)/y"'-; then gi{y) strictly decreases, hi{y) strictly in- 
creases, and the the fraction g'{y)/h'{y) — gi{y)/hi{y) strictly decreases as 
required.) Then, in particular, the function r + g{y) — a^h{y) is unimodal, and 
there exist the single positive root i] of the equation r + g{y) — a*/i(y) = 0. 
Therefore the function ([57)1 is unimodal for < y < 77. The branch x = X{y), 
< y < y*, is repulsive, the branch x — X, y^ < y < i], is attractive, and 
(a;,,y,) is the unique turning point. 

For a ~ a, the system may have two types of canards, see Figure [T51 Early 
canards may exist for xo,yo > satisfying ^(yo) < xq < ^(0), and they have 
the standard structure. Late canards, which may exist for a;o,yo > satisfying 
< yo < y*: Xf, < xo < X(yo), are more interesting. A late (xq, yo)-canard 
exhibits additional delayed loss of stability phenomenon: after following down 
at a; « xq it follows closely the axis y = until a point x !v S^q > q/r, and then 
jumps up to the attractive branch of the slow manifold. The point ^0 is the 
solution of the equation ^'^exp(— r^) — a;Q exp(— ra;o). Indeed, close to the axis 
y — the dynamics is governed by the equation dy/dx = y{q — rx)/{px) whose 
solutions satisfy the relationship Iny^' — In a;' -|- ra; = const. 
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Figure 13: Early and late canards for the modified Lotka-Vilterra system. Early 
canards may exist for a;o,yo > satisfying X(yo) < a^o < -'^(O), and they have 
the standard structure. Late canards, which may exist for a;o,?/o > satisfying 
< yo < y*: < xq < X(yo), are more interesting. A late (iq, yo)-canard 
exhibits additional delayed loss of stability phenomenon: after following down 
&i X Ki xq it follows closely the axis y = until the point and then jumps up 
to the attractive branch of the slow manifold. 

Proposition 6.5. There exists an early (xq, yo) -periodic canard for any xq, yo > 

satisfying X(i/o) < xq < ^(0), and there exists a late {xo,yo) -periodic canard 
for any xo,yo > satisfying < yo < y*, x^ < xq < X{yo). 

Proof. As a± we choose any numbers which are sufficiently close to a, and 
satisfy a_ < a* < a+. Note that 



and the inequalities (|5D|) and ([5T|) follow. The construction of the region D in 
the case of an early canard is the same as in the first example, and in the case 
of the late canard is explained in Figure [T31 Here ^» denotes the unique positive 
solution of the equation ^'?exp(— rf) = a;* exp(— rx*). Nonexistence of confined 
in D [J r cycles for a — a-, a+ can be proven as in the previous examples. 
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X 



Figure 14: The region D bounded by the curve F. At last one cndpoint of the 
velocity vector points strictly outward of D at all points of F, except of two 
points: one is our "target point" {xo,yo), and another is the bold point at the 
eastern bound of D. However, there are no cycles which touch F at the second 
point, since the longest possible travel along the axis y = ends at the point 
2q/r — x*, which is located strictly to left of the eastern bound of D. Thus, by 
Proposition 16. 1[ there exists a cycle which is confined in DljF, and touches F 
at {xo,yo). 

Proof of Proposition 16.11 

The Poincare-Bendixson theorem can be stated in several ways. The statement 
that is relevant to the equation (j49p is the following. Suppose S is a closed, 
bounded subset of the plane; S does not contain any fixed points; and there 
exists a trajectory confined in S. Then either this trajectory is a closed orbit, 
or it spirals toward a closed orbit. 

For a particular value of a a solution x{t) of ([15)) is called directed, if x(0) G F 
and x(i) G D for i > 0. There exists a directed solution x(t) for a = a+,a_. 
To prove this claim we suppose that tr Ja_ < 0, and consider a solution which 
begins in a sufficiently small vicinity of ea_. Then \y{t) — ea_ | — s- as i — > oo, 
and \y{t)\, is bounded from below at t < (because e^^ is a sink due to 
det Je„ jtr < 0). By the Poincare-Bendixson theorem y{t) must leave D 
in negative time (because there is no cycles at a = O-); in particular, y{t) 
touches F for the first time at some t = t < 0. It remains to set x(t) = y(f + r). 
Analogously, using the backward time, we prove that there are no T-directed 
solutions at a — a^. 
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Denote by ao G [a_,a+) the upper bound of a € [a^,a+) for which there 
exist some directed solutions. For a — gq a directed solution Xo(t) also exists 
by limit transition. If Xo(-) is periodic, then the proposition holds. To finalize 
the proof we suppose that x(-) is not periodic, and arrive at contradiction. 

By the Poincare-Bendixson theorem there are only two possibilities: either 
(a) Xo(f) spirales toward C a cycle C C f, or (b) Xo(t„) — > for some tn — > oo. 




Figure 15: The trapping curve Fq. 

Let Fq C -D be a Jordanian curve which bounds the open domain Dq, and r 
be a positive number. We say that the pair {Fq, r} is trapping if simultaneously: 
the set Dq [J Fg is forward invariant for the equation x — f(x, ao), and x(r) G Dq 
holds for any solution satisfying x(0) G Fg. If a trapping pair exists, then for a 
slightly greater than ag the solutions of equation (|^ that begins at Xo(0) are 
also attracted to arbitrary small vicinity of Dq. That is, there exist a directed 
solutions at some a > oq. Thus, to arrive at contradiction it is enough to 
construct a trapping pair. To this end in the case (a) we choose a point y G C 
and consider the corresponding outward normal n to C. Let A satisfies the 
relationships [y,y + An] C D, and f(y + An, ag) • f(y, ag) > 0, < A < A. By 
definition, the solution xo(i) crosses the segment (y,y + An] infinitely many 
times, see Figure [15] Let to and ti be two successive moments of such crossings 
with the corresponding values Ag, A2. Consider the curve Fg which consists of the 
trajectory Xg(t), tg < i < ti, together with the segment [xo(io), Xo(ii)]. Since 
Xo(t) spirales toward C, the inequality Ag > Ai holds. Therefore, {Fg,ti— <g + l} 
is a trapping pair, and we arrived at contradiction in the case (a). 

By tr J{eag ) < the case (b) can be partitioned in turn into three cases: (bl) 
is a source; (b2) is a sink; (b3) e^g is a center in the linear approximation. 
In the case (bl) we immediately arrive at contradiction with the condition (b). 
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In the case (b2), for a slightly greater than uq, the solution which begins at Sag 
is attracted to a small vicinity of eQ^. Thus, there exist directed solutions for 
some a > ao, which contradicts the definition of oq. It remains to consider the 
case (b3), which is similar to the case (a) above. Indeed, consider a segment a = 
{eao,eao + z] where z is close enough to Bag to guarantee that a C D, and that 
f (Yi flo)j y G cr, is not coUinear to z. (This can be done because Bag is a center 
in the linear approximation.) By the condition (b) the solution Xo(i) crosses 
the segment a infinitely many times. Let to and ti be to successive moments 
of such crossings. Consider the curve Tq which consists of the trajectory Xo(t), 
to < t < ti, together with the segment [xo(to), Xo(ti)]. By construction the pair 
{ro,ti — to + 1} is a trapping pair, and we arrived at contradiction in the case 
(b3). The proposition is proven. 

7 Non-smooth perturbations 

Consider a perturbed system ([5^ : 

i^y, £y = -X + F{y - a) + F{x,y,a), 

where F is continuous and small in the uniform norm: sup \F{x, y,a)\ < 5 <^ 1, 
but there is no bounds for on its derivatives. In this case applicability of usual 
tools is doubtful. 

Proposition 7.1. There exist e,6 > such that for < e < £,0 < 5 < 5 there 
exists a small and a periodic canard{x^{t)^y^{t)) of the system x = y, sy = 
—X + F(%j — Ce) + F(x, y, ae) which satisfies m.aKt{xe{t)} = b. The trajectory of 
this canard approaches T{b) as e^5 —>■ 0. 

Proof follows from the following modification of Proposition 16.11 Consider the 
equation 

X = f(x,a) +f(x,a). (68) 

Here f(x, a) is continuous and uniformly small: supf(x, a) < S ^ 1, but there 
is no restriction on its derivative. Under conditions of Provosition [KT\ for some 
a £ (a_,a+) there exists a cycle of system (j68p which belongs to D[JT, and 
which touches T. Proof is essentially the same as of Proposition [^Tj 
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